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Abstract 

We study a free quantum motion on periodically structured manifolds 
composed of elementary two-dimensional "cells" connected either by 
linear segments or through points where the two cells touch. The 
general theory is illustrated with numerous examples in which the 
elementary components are spherical surfaces arranged into chains in 
a straight or zigzag way, or two-dimensional square-lattice "carpets" . 
We show that the spectra of such systems have an infinite number of 
gaps and that the latter dominate the spectrum at high energies. 



1 Introduction 

The spectral behaviour of periodic systems is of a great importance. Having 
typically a band structure, such spectra differ by the number and structure 
of the gaps. For usual Schrodinger operators the number of gaps is generi- 
cally infinite in the one- dimensional situation and finite in higher dimensions. 
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Moreover, the gap widths decrease as the energy tends to infinity, the rate 
of decay being tied to the regularity of the potential. 

In case of a singular periodic interaction the gaps may not close. A canon- 
ical example is the Kronig-Penney (KP) model, i.e. a chain of 5-potentials 
where the gaps are asymptotically constant |[AGH|1 . Even more singular cou- 



plings like generalized point interactions may exhibit gaps which are growing 
at the same rate as the bands [ EGrj| , or even grow while the band widths are 
asymptotically constant. A typical example of such a behaviour is a modifi- 
cation of the KP model with a chain of the so-called ^'-interactions IIAGH 



This behaviour is not restricted to one dimension; similar results can be de- 
rived e.g. for lattice graphs with appropriate boundary conditions coupling 

For instance, the cor- 



the wave functions at the vertices 



Large gaps has interesting physical consequences, 
responding Wannier-Stark problem in which we add a linear background 
potential to a periodic chain of ^'-interactions has counterintuitive proper- 
ties: the absolutely continuous spectrum of the corresponding Hamiltonian 
is empty ||AEL|| , and in fact, the spectrum is known to be pure point for 
"most values" of the potential slope ||ADE|| . These results can be explained 
by observing that tilted gaps represent classically forbidden regions and that 
their large widths prevent the particle of propagating over long distances. 

On the other hand, the physical meaning of the (^'-coupling remained un- 
clear for a long time. Recently it has been demonstrated that this interaction 
can be approximated in the norm resolvent sense by a family of Schrodinger 
operators - see ||ENZ|| where also a bibliography to the problem is given - but 



previous studies brought some interesting non-potential approximations. An 
interesting example is given by a "bubble scattering" in which two halflines 
are attached to the surface of a sphere - see |[Ki|| and also |ETV|, Br|. Such a 



system typically exhibits numerous resonances but the background transmis- 
sion probability dominates and vanishes at the limit of large energies. This 
observation is of importance because systems of a mixed dimensionality are 
not just a mathematicians toy, but they can model real objects such as a 
fullerene molecule coupled to a pair of nanotubes | [Ka| .[] 

The aim of the present paper is to study systems with components of 
different dimension in the periodic setting. We intend to demonstrate that 
the structure of the configuration space in this case is reflected in the gap 



1 Another model for such systems could be that of manifolds connected smoothly by 



thin tubes. Existence of gaps in this setting was demonstrated recently by Post [Po 
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behaviour. After describing a general method to couple periodic systems 
of spheres, either joining them by line segments or directly through points 
where they touch, we will discuss in Sections 3-5 a number of examples. 
The results, summarized in Proposition |5.1| , show that in all the considered 
cases the number of gaps are infinite and the gap-to-band width ratio in- 
creases with the band index. The estimated growth is slower than in the 
case of the ^'-interaction, and it is slower for a two-dimensional lattice than 
for a linear chain, but it is still powerlike for spheres joined by linear seg- 
ments, thus confirming our conjecture that the effect is related to the change 
in dimensionality the particle must undergo. Even for a tighter coupling, 
however, where the spheres are coupled directly through contact points, the 
gap-to-band ratio is still logarithmically increasing. 



2 General theory 

2.1 Building blocks of the Hamiltonian 

Suppose that Xq is a two-dimensional Riemann manifold. By Hq we denote 
a Schrodinger operator, 

#o = W^i-idj - A^gl^^i-idk -A h )+V, 

on L 2 (X , Ig^^dx) with smooth vector and scalar potentials. The formalism 
we are going to describe extends easily to the case dimX = 3 but we will 
limit ourselves here to referring to ||BG2|| for guidelines concerning such a 



generalization. The metric structure of Xq is fixed and we will employ the 
shorthand notation L 2 (X ) for simplicity in the following. Let further Xj,j = 
1, . . . , n, be a finite or semiinfinite line segment which can be identified with 
the interval [0,dj), < dj < oo. No external potentials are supposed to act 
on the particle on Xj, i.e. we consider the free operators Hj = — d 2 /dx 2 
on L 2 (Xj) with Neumann's condition at the endpoints (the "right" endpoint 
x = dj requires a boundary condition only if dj < oo) as the building blocks 
of the system Hamiltonian. 

As we have said above we consider systems with configuration space con- 
sisting of infinite number of copies of a manifold which are connected either 
by isolated points common for the pair of neighbouring copies, or by line 
segments connecting such points. We will concentrate on the latter case 
which is more complicated. The former one can be regarded as the limiting 
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situation where the length of the connecting segments tends to zero, and the 
corresponding modification of the formalism is easy. 

A building block of our model is thus a "hedgehog manifold" obtained 
by attaching each segment Xj to X at a point qj G X , or more exactly, by 
identifying the point G Xj with qj G X ; we suppose that all the connection 
points qj are mutually different. The topological space constructed in this 
way will be denoted as X; it can be endowed with a natural measure which 
restricts to the Riemannian measure on Xq and to the Lebesgue measure on 
each Xj , j — 1, . . . , n. This yields the identification 

L 2 (X) = L 2 (X ) © L 2 (Xi) © • • • L\X n ) 

for the Hilbert state space of the system. 

By So we denote the restriction of the operator H defined above to the 
family of functions 

{feV(H ) : /( gi ) = --- = /(<7n) = 0}, 

which obviously makes sense as long as dimX < 3. In a similar way 
we use the symbol Sj, j = l,...,n, for the restriction of Hj to the set 
{/ G T>(Hj) : /(0) = 0}. The Schrodinger operators we consider are by defi- 
nition self-adjoint extensions of the symmetric operator S = So ©Si©- • -ffiS n ; 
their construction is a standard matter discussed in numerous papers starting 
with |[ES1|| . The most efficient way to describe them is based on a bijective 
correspondence with the Lagrangian planes in Q x Q, where Q = C 2n . To 
describe it, we introduce the boundary- value operators 



r 1 , T 2 : V(S* 



by 



T\f) := (a(/o,gi),...,a(/o,g n ),-/i(0),...,-^(0)), 

T 2 (/) := Wo, ft), ... , 6(/o, Qn), /i(0), • • • , /„(0)) . (2.1) 

Here a(f ,qj) =: aj(f ) and b(f ,qj) =: bj(f ) are the leading-term coeffi- 
cients of the asymptotics of fo in the vicinity of the point qj as determined 
|BG1| , |BG2|| , or the generalized boundary values 27r(dimX — 1)L and Li, 



in 



respectively, in the terminology of |[ES1] , [ES2| 



Let A be a Lagrangian plane in Q x Q, i.e. A 1 - = A with respect to 
the skew-Hermitean product [x\y] := (£i|2/2) — (-^2 1 2/i ) in Q x Q. Then any 
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restriction of the adjoint operator S* to a family of functions from T>(S*) 
specified by the boundary condition (r 1 /, T 2 f) 6 A is a self-adjoint operator 
which we denote H A . Recall that a Lagrangian plane is, in general, the graph 
of a self-adjoint operator L : Q — ► Q so that the above boundary condition 
can be rewritten as T 2 f = L(r 1 /). To avoid problems with the invertibility 
of L one can view A also as the graph of a "multivalued" operator in Q, in 
other words, one may describe it through a relation Lx = My, (x, y) G Q xQ, 
where L,M : Q — > Q are linear operators satisfying the conditions ||KS|| : 

(i) LM* = ML*, 

(ii) rank(L, — M) = n . 



2.2 The resolvent 



We are concerned with spectral properties of the said self-adjoint extensions, 
which are as usual defined from the resolvent. The latter is expressed here by 
Krein's formula ||AGH| , App. A]: If we denote by H° the decoupled operator 
Hq © Hi © • • • H n , then we have 



{H 1 



7(*)[Q(*)-A]-Y(s) 



(2.2) 



for any z in the resolvent set, in particular for z ^ C \M, where the operator 
l{z) : Q — > Ti is given by the formula 



70) := (r 1 \N Z ) \ N z — Ker (S* —z) 
and Q(z) : Q — > Q is defined as 

Q(z) := T 2 7 (z) . 

Then the inverse [Q(z) — A] -1 exists for all non-real z. To find an explicit 
expression for the Green function of the operator H A from ( |2.2|) we need to 
know the Green function Go of H°. 

Notice first that it is easy to find the Green function Gj of Hj: one has 



Gj(x, x'; z) 



cosh [—y/—z(dj 



'])]+ cosh [-y/=z(dj - (x + x')] 



2v r —z sinh [— yj— zdA 



(2.3) 
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Using the natural decomposition Q = C 2n = C" x C" we write the matrix 
representation of the operator [Q(z) — A] -1 in block form, 



[Q(z)-A]- 1 



T{z) W(z) 
U(z) V(z) 



(2.4) 



Since H = Ho © Hi © • • • H n) the Green function of the operator H A can 
be represented as a matrix of integral kernels of operators acting from Hk to 
Hj, 

G A (x, x'- z) = {G%( Xj , x' k - z)) ^ n with Xj E X 3 , x' k E X k . (2.5) 
Let (£1, . . . , £„, 771, . . . , rj n ) G Q = C n x C n , then a direct calculation yields 

■ ■ ■ ,€n,Vl, ■ ■ ■ ,Vn) 
' n 

This implies the adjoint operator action, 

7*(^)(/o,/l, •••,/n) = (tl,---,tn,Vl,---,Vn) , 

where 

& = / Go(qj,x;z)fo(x)\g(x)\ 1/2 dx, r}j = Gj(0,x; z)fj(x) dx . 

J Xq J Xj 

The matrix Q(z) then has block-diagonal form 

Q(z)~- 



Q ll {z) 
Q 22 (z) 



here Q u (z) coincides with the Q-matrix Qq for the pair (Sq,Hq). Recall 
that the Q-function in the Krein formula always corresponds to a pair of a 
self-adjoint operators and fixed symmetric restriction. In the present case it 
has the form 

Qi\z)=GT\q v q k] z), (2.6) 

where G r Q en is the renormalized Green's function obtained from Go by sub- 
tracting the diagonal singularity, 

G (x , x' ; z) , ifx ^x' ; 
Gr(x ,x' ;z) = { Um 



G (x , y ; z) + — In p(x , y ) 
Ztt 



, if x = x' . 



Here p(x ,y ) denotes the geodesic distance on X . On the other hand, the 
matrix Q 22 (z) is diagonal, 

Qf k (z) = 6^0,(0,0; z). 

Using the above formulae we can write the matrix element kernels in (|2.5|) 
more explicitly, 

Gjk(%ji ^jkGj^Xji Xj, z) K j k i%j i % k i i 

where 

n 

K 00 (x , x' ; z) = ^ t jk (z)G (x , qf, z)G (q k , x' ; z) , 
j,k=i 

n 

K ok (x , x' k ; z) = G k (0,x' k ;z)^2w jk (z)G (x ,q j ;z) , k> , 

i=i 

n 

K j0 (xj,x' ;z) = G j (x j ,0,;z)^2u jk (z)G (q k ,x' ;z) , j>0, 

k=l 

K jk (xj, x' k ; z) = Vj k (z)Gj{xj, 0, ; z)G k (0, x' k ; z) , j,k>0; 

the coefficients refering to the block representation (|2.4j ), (tj k (z)) = T(z), 
etc., can be in principle computed explicitly. 

2.3 Coupling hedgehog manifolds 

In the next step we are going to glue together the building blocks considered 
so far. To begin with, we consider such a manifold X and select some number 
of finite segments of lengths d±, . . . , d s , 1 < s < n. At the same time, we fix 
a finite number of distinct points pi, . . . ,p m G X such that {p±, . . . ,p n } fl 
{qi, . . . , q n } = 0. We fix a Hamiltonian H A on X and consider its restriction 
S to the set of functions 

{ / G V(H A ) : f( Pl ) = ■■■ = f(p m ) = fid,) = ■■■ = f{d.) = } . 

Let us find the Q-matrix of the pair (S*, H A ) which is a (m + s) x (m + s) 
matrix Qiz) with block structure, 

Q l \z) Q 12 iz)- 
Q 21 {z) Q 2 \z) _ ■ 
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Using the formula for the Green function of H A we can write the elements 
of the above matrix as 



Q 



Q 



S21 



Q 



-Koo(Pj,p k ;z) , 
-K ok (pj,d k ;z) , 

-K j0 (dj,p k ; z) , 



1 - S jk )Go(pj,p k ;z) 
1 < j, k < m , 
1 < j < m , 1 < k < s 

1<j<s, l</c<m. 



(2.7) 



22 1 



SjkGjidj, dj\ z) - K jk (dj, d k ;z), l<j,k<s. 



Recall that Gq 611 denotes the renormalized Green's function; we drop of course 
the superscript whenever the two arguments are different. 

The coupling will be realized through conditions relating the generalized 
boundary values. We will not strive for at most generality, however, because 
formulae encompassing manifolds with arbitrary n, m would be rather cum- 
bersome. We will instead discuss in some detail properties of a quantum 
particle living on chained manifolds of different dimensions, i.e. the case 
m = n = 1; later on we will extend the argument to a particular situation 
with m = n = 2. 

Consider, therefore, a manifold Xo on which a pair of mutually different 
points p, q is selected. At q, a segment of a length d is attached, while p is 
a "socket" to which another "tailed" manifold can be coupled. In analogy 
with (12.61) we introduce the matrix 



Qo(z) 



G^(q,q;z) G (p,q;z) 
G (q,p;z) G v en (p,p;z) 



(2.8) 



and similarly, the segment will be characterized by 

Qi(z) = 



Gi(0,0;z) Gi(0,d;*) 
Gi(d, 0;z) Gi(d,d;z) 



(2.9) 



Using (|2.3| ) we find 



6* 



or 



coth (\/—zd\ + 



Ql^ z ) = 6 ^ cot (kd) + 



1 - s jh 



^z sinh (\/— zd) 
1 - Sjk 



k sin(kd) 



(2.10) 
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in the usual momentum notation, k := i\J —z for z£C \ R + . 

The operator H A on X is specified by the boundary conditions at the 
point q identified with the left endpoint of the segment, £ [0, d). In general, 
these conditions can be given in the form, 



b(fo,q) = <xf{(p) + Pa(fo,q) 
A(0) = lfi(0)-aa(f ,q) 



(2.11) 



with /3, 7 £ R and a £ C; we suppose a 7^ such that the manifold X and 
the segment are coupled in a nontrivial way. For the sake of simplicity we 
will restrict ourselves to the case where f3 = 7 = 0, i.e. 



Kf ,q) = afi(0), /i(0) = -aa(/ o ,g) 



(2-12) 



This can be regarded as a "minimal" coupling between the two configuration- 
space components, because in the "switched-off state", a = 0, the manifold 
Hamiltonian contains no point interaction at the point q and the segment 
part satisfies the Dirichlet condition at x% = 0. Notice, however, that there 
are other natural choices such as 



a 



7T 



-7r(l + m^), 7 = 2p, 



which describes the particle passing through the junction at a low energy if 
the segment models a thin tube of radius p - cf. [|fcCS2| . 

The boundary condition ( 2.12|) can be cast into the form given in Sec. |2Tl 
if we choose M as the 2x2 unit matrix and 



L :-- 



a 
a 



(2.13) 



The Q-matrix entering Krein's formula for the operator H A can be expressed 
in terms of the matrices (12.81) and (t^.91) as 



Q(z) 

From (|2~T3|) and ( [Hp we find 

[Q(z) - A]- 1 - 



Ql 1 





Ql 1 



(2.14) 



Ql\z)Q?{z)-\a\ 



Q\\z) -a 
-a Q l l (z 
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and therefore 



G 00 {x , x' ; z) — G Q (x Q ,x' ; z) — u >„ n/ 7 — ; — — Go(xo,q; z)G Q (q,x' ; z), 



Q${z)Q?{z)-\a 



a 



G m (x , x[; z) — n ,~ u , 7 — ; — rr7 G (x 0} q; z)Gi(0, x[; z) , 

Qo {z)Qi_\ z )-\ a \ 

G^ (x 1 ,x' ; z) = „ u , , nll , > — —^G 1 (x 1 ,0;z)G (q,x' ;z), 

Gii(xi, x'i, z) = G\(xi, x'i, z) — n \niu \ — i — 12 ^i(* E i' 0> z )Gi(®, x[; z) . 

Thus we can calculate the matrix elements fl2.7| ) (the indices j, A; are trivial 
in the present example and we will drop them): 



2 



^ *ggg^ , (2-15) 



g^g^-lal 2 

gg^)gf(z) 
g^g^Hal 



S22^ x Q^QPWfO*) 



G"{z) = Qf{z)- 



Q£{z)Q\\z)-\a\* ■ 
These formulae can be made even more explicit by plugging in ( p,10| ). 



2.4 Point-coupled manifolds 

In the same way one can treat the limiting situation when the lengths of the 
connecting segment shrink to zero. Then only the boundary conditions have 
be modified. Consider the simplest case when X and X 1 are coupled by 
identifying the points pj G Xj, j = 0,1. The generalized boundary values 
(|2.1| ) are then replaced by 

r^/oj/i) := Wo,po),a(fx,Px)) , 
r 2 (/o,/i) := (6(/o,Pb),6(/i.Pi)) • 



Such a coupling was first discussed in | |ES3| ] in the situation where X and X 



-1 



are two planes. The four-parameter set of all possible self-adjoint extensions 
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was described there and the result adapts easily to more general manifolds. 
For the sake of simplicity, however, we will again restrict our attention to the 
"minimal" coupling given by the conditions 

Kfo,Po) = ota(fi,pi) , 6(/i,pi) = aa(f ,po) (2.16) 
with a complex parameter a, decoupled manifolds corresponding to a = 0. 

3 Infinite necklaces 
3.1 General periodic case 

As an illustration of how to couple "hedgehog" manifolds, we are now going 
to analyze now the simplest nontrivial case i.e. when the building blocks dis- 
cussed above are chained into an infinite "necklace" . To define the Hamilto- 
nian we have to specify the boundary conditions coupling the outer endpoint 
of the segment of the first building block, starting at q, to the point p of the 
second one. The boundary-value operators T 1 and T 2 for the operator S are 
of the form 

^(/o,/!) := (a(f Q ,p),f[(d)) , 
r 2 (/o,/i) := (b(f ,p),h(d)) . 

Notice the positive sign of f[(d) in comparison with (|2.1|) which reflects the 
orientation of the segment [0, of]. 

Consider now a countable family of identical copies of the manifold X, 
i.e. Xm = X for all m 6 Z and set Z := {J me %X m . The state Hilbert space 
of this necklace is 

oo 

L\Z)= L\X m ). 

m=—oo 

Schrodinger operators on the necklace will be identified with self-adjoint ex- 
tensions of the symmetric operator S := © meZ S m , where S m := S for any 
m G Z. Obviously, the boundary-value space of S is of the form 

oo 

Q = Qm with Q m = C 2 forall m 
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and 



f j = Y j m with f 3 m = f j forall m and 3 = 1,2 



m=—oo 



Of course, the operator S has infinite deficiency indices, and therefore plenty 
of self-adjoint extensions. We restrict our attention to those which are local 
in the sense that exactly the point d of X m is coupled with the point p 
of X m+ i. Moreover, we will consider the situation when the coupling d to 
p and q to is minimal in the sense described above. Consequently, for an 
element g = {g m } G Q with g m = (fo >m , /i, m ) we impose boundary conditions 
analogous to Q2.12|) : 

Kfo,m,p) = a/v m _i(0) , /i,m(0) = -ao(/ 0>m+ i,p) ; 



this can be written concisely as 



t 2 g = Lt l g, 



(3.1) 



where L is an operator in Q given by a matrix L = (L mn ) m>ne z, where 
L mn = if \m— n| 7^ 1 and 



a 




J rn+l,rn 





a 



We then infer that the self-adjoint operator H L specified by the boundary 
conditions ( |3.1| ) has the following resolvent 

(H L - z)- 1 = (H° - z)- 1 - i{z)[Q{z) - L^Tiz) , 

where Q(z) = {S mn Q(z)}. In this way, the dispersion relation for H L can be 
obtained by introducing the quasimomentum 9 G [0, 2tt) and performing the 
Fourier transformation of the operator Q(z) — L. Thus result is an operator 
in the space L 2 ((0, 2w)) <S> Q with kernel 



P(9,z):= J2 (Qmo(z) - L m0 ) e me = Q(z) - 

m=— 00 

The dispersion relation is of the form det P(9, z) = 0, or 
det 



a e 



a e 
w 



w 



Qn(z) 

Q 2 i(z) - ae~ 



Qi 2 (z) - ae" 
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which is equivalent to 

detQ(z) - (Q 12 (z)ae- te + Q 21 (z)ae ie ^ - \a\ 2 = 0. (3.2) 

As in similar situations, we have isospectrality with respect to the coupling- 
constant phase: put tp = arga, i.e. a = \a\e lip , then the last condition can 
be written as 

det Q(z) - \a\ (Qi 2 (z)e-< e+ ^ + Q 21 (z) e^ + ^) - \a\ 2 = 0, 

which shows that without loss of generality we may restrict ourselves to 
the case a > 0; this we shall assume in the following. Using the fact that 
Q 2 i(z) = Qu(z) holds for real z, the condition ( |3.2|) can be rewritten as 

detQO) - \a\ 2 = 2a (ReQ 12 (z) cos 8 + lm Q i 2 (z) sin . (3.3) 

Hence a necessary condition for z G spec (if) is 

| det Q(z) — \a\ 2] 



2a\Q 12 (z) 



< 1 • (3.4) 



If Q\ 2 (z) = Q 2 \(z), which is true in particular if H° is a real operator (i.e. 
commutes with the complex conjugation), the relation ( |3.3j ) simplifies to 

det Q(z) - \a\ 2 

cosO = , 3.5 

2aQ l2 {z) 

and the condition (|3.4|) becomes necessary and sufficient. If H° is real, the 
condition ( |3.5|) can be made more explicit: using (|2.15| ) and the fact that 



Qj = Q 2 holds for j = 1, 2, we find after a short computation 

det Qq(z) det Q 1 (z) - 2a 2 Q 1 1 (z)Q\ 1 (z) + a 4 
~ 2a*Ql 2 (z)Q\ 2 (z) 

Furthermore, using ( |2.10| ) we get 



det Q (k 2 ) sm(kd) - 2a 2 k cos(kd) Ql\k 2 ) - a 4 k 2 sin(kd) 

2a 2 kQ l2 {k 2 



UCO mO\rh J OLLL^IXU, J ZiVL IX UUO^/vU J l^Q {IX J LX IX SLLL^IXU, J . 
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3.2 Spherical beads 



Since our aim is to present solvable examples, we study next the situation 
when the elementary building-block manifold Xq is a two-dimensional sphere 
§ 2 of a fixed radius a > 0. We parametrize it by spherical coordinates, 



x 

y 

z 



a cos i? cos if . 
a cos i? sin ip , 
a sin $ , 



with d G [— vr/2, vr/2], ip G [0, 27r). We will assume that there are no external 
fields, so the starting operator for construction of the Hamiltonian is the 
Laplace-Beltrami operator Alb on S 2 . Its Green's function is an integral 
operator with the kernel 



G (x,y; z) 



4 cos(7rt) 



— cos 



p{x,y) 



(3.7) 



where V\ is the Legendre function, p(x, y) is the geodetic distance on the 
sphere, and 

t = t(z) := ^Vl +Aa 2 z. 
This allows us to express the renormalized Green's function, i.e. we find 



l 

2^ 



7T 



tan(7rt) — In 2a + Ce 



(3.8) 



(see e.g. [ [BE , Tab. 3.9.2]), where Ce is Euler's number and ip the digamma 
function. We use again the conventional notation z = k 2 for the energy 
parameter; if there is no danger of misunderstanding we will often supposes 
the dependence of various quantities on k. 



3.3 Loose necklaces 

We shall next consider two particular segment-connected periodic chains: 

Example I: Suppose that the connecting segments are attached at antipodal 
points as sketched in Fig. 1 so that the geodesic distance of the junctions is 
7ra. We will denote the segment connecting the spheres § 2 and S 2 +1 as J n , 

with the endpoints 0^ = pf ] G S 2 and = p^ +1) G § 2 +1 . The lower- 
index numeration is somewhat arbitrary and serves just to having a common 
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Figure 1: A loose straight necklace 




Figure 2: A loose zigzag necklace 



notation for the present configuration and that considered below. 

Example II: Alternatively, assume that the junction points are chosen on one 
pole and on the equator point, as sketched in Fig. 2, such that their geodesic 
distance is na/2. The segment I n now connects the points = p± £ §^ 
anddW=^ n+1) G 



While the diagonal part ( |3.8| ) of the matrix Q does not depend on the way 
we arrange the spheres, the off-diagonal parts differ and now become 



^0 
^0 



"k r ( 



f + |) cosvr(I + |) 



4 cos(7ri) 



(3.9) 
(3.10) 
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K<2 R<2 r;2 

Figure 3: A tight straight necklace 

for the zigzag and straight case, respectively, with the notation we have 
adopted. In the same way, the dispersion relation (|3.6|) becomes 

o lj 

Ql l Ql l ~ QoQo - 2a 2 k cot(kd) Ql l - 2a 2 k . Vo cos 9 - a 4 k 2 = (3.11) 

sm( Kd ) 

with j = 2, 3 in Examples II and I, respectively. Let us remark that the 
condition with j = 2 is valid whenever all the Q^ 2 are the same. Hence the 
spectrum does not change when we rotate an arbitrary semi-infinite part of 
the chain around the axis given by the appropriate connecting segment, such 
that, geometrically speaking, the zigzag chain need not be periodic. 

3.4 Tight necklaces 

In a similar way, one can treat periodic sphere chains which are connected 
through points where they touch (i.e. shrinking the line segments to zero), 
with the boundary conditions (|3.1| ) replaced by (|2.16|) at each junction. We 
shall consider again two particular situations analogous to the periodic chains 
discussed above: 

Example III: Suppose that the junctions are situated at antipodal points as 
sketched in Fig. 3, being obtained by identifying the points p± G and 

(n+1) -2 
P3 fc s n+l- 

Example IV: The tight zigzag chain in Fig. 4 is obtained by identifying the 
points pi n) G §^ an d p ^ e S£ +1 . 
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Figure 4: A tight zigzag necklace 



The dispersion relation now reads 



Ql'Ql 1 - QoQo - 2« Ql j cos e + a 2 = 



(3.12) 



with j — 2,3 corresponding to the Examples IV and III, respectively. 



So far we have considered only "manifolds" with a linear structure. Having 
in mind essential differences between spectra of periodic Schrodinger oper- 
ators in different dimensions to detect, it is useful to also to look systems 
which are periodic in more than one direction; we will do this again by first 
analyzing simple examples. This time we arrange our spherical "beads" into 
a square lattice, coupling them either by line segments or directly through 
touching points. 

Example V: Suppose that the connecting segments are attached at four 
equally spaced points at the sphere equator as sketched in Fig. 5, where 
the labeling of the junctions and segments is indicated. With the notation 
introduced in Fig.5 the boundary conditions defining the Hamiltonian read 



4 Square bead carpets 




= —a 



= aa 



= —a 
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Figure 5: A loose square bead carpet 



b(ft m \p { 2' m) 



= aa 



( f(«>™) 

yo 1P3 

-a [fl ) 

aa (ft' m \p ( 2' 
-a (yf™-*)) 
aa ( / ( 



(n,m) (n,m 



g(n,m+i)j 
g(n,m-i)j 



The dispersion relation is derived as in the previous section, but it becomes 
rather cumbersome. It is useful to introduce the following notation: 

A := 1 (Q W - QW) + (Qo 1 oos(Jbd) + Q 2 ) + a 4 , 



6i 



+ 2 Qo Qo 



k 2 sin 2 (fed) 
1 



ij+i 



k 2 

Ql' 2 ~ J 



cr 



fcsin(fcd) 
a 2 



/c 2 fcsin(A;d) 



, j = o, i , 



cos(fcd) + a 2 Q 1 1 



k 2 sin 2 : (fcd) 



fc sin(fcd) 

Qo^Qq 1 ~~ Qo 2 Qo 2 
k sin(fcd) 



A cos(fcd) 
A; sin (fed) ' 



fc sin(/cd) 



2 n l,4-j , Qq' 4 ^Qq 1 cos (kd) 



fc sin(A;d) 
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Figure 6: A tight square bead carpet 




i = i,2. 



Using this notation, we can write the spectral condition as 



(a 2 - a\){bl - b\) + {c\ - c 2 ) 2 - 2[(ci + c 2 ) 2 (a b + a^) 
-2cic 2 (a + ai)(6 + h)] 

+2aA [(a b + ai&i - c\)ci + a bi + ai& - cic 2 )c 2 ] (cos 9 1 + cos 2 ) 
+2a 2 A 2 [(c 2 - c 2 ,) cos(^i + 2 ) + (c 2 - ai&i) cos(0i - 2 ) + c 2 - a &o] 



where #i,0 2 are the quasimomentum components. 

Example VI: This arises from Example V by shrinking the connecting seg- 
ments to zero, as indicated in Fig. 6 where the labeling of the junctions is 
the same as in the previous example. After a straightforward calculation we 
find that the spectral condition now takes the form 



2a 3 A 3 Ci(cos#i + cos# 2 ) + a 4 A 4 = . 



(4.1) 




(4.2) 
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5 Gap dominance at large energies 



As customary in periodic systems the spectrum of the above described opera- 
tors (which we denote by Hi, ... , Hyi according to the example number) has 
band structure. To see how the gap width and the band width are related at 
high energies, consider first the points 



k n :— — , k" n :— — , n — 1, 2, ... , (5-1) 

for which sm(dk' n ) = cos(dk' n ) = cos7r (| + §£(&2n)) = 0' sucn that the func- 
tions <5q and Qi have poles. Thus it is natural to look for spectral bands in 
the vicinity of these points. We fix e > and denote by J' n = [k' n — 5' n , k' n + S' n ] 
the maximal closed neighbourhood of the point k' n in which the inequality 

| sm(kd)\ < k~ e 

is satisfied; in the same way the intervals J" = [k" — 5",k" + 5"} and J"' = 
Win ~ ^2n + K'] correspond to the inequalities 



cos(kd)\ < k e and 



(\ 




COS 7T h 





\4 





<k~ e , (5.2) 



respectively. It is clear that all the 8' n , . . . ,5'" are strictly positive, and it is 
not difficult to check that 

S'n ~ d-\k' n r , C C ~ 2(na)-\k' n r , c C ~ ^a)-\k' 2n r 

as n — > oo. Our aim is to show that for a sufficiently high energy the spectral 
gaps contain the complement of the above intervals. More specifically, define 

n K := [K,oo)\{J(j> n uji;uj:) 

71=1 

for a fixed K > O.In this set, we have (sm(kd))^ 1 = 0(k £ ) as k — > oo, and 
similarly 

Qfc=0{lf), Ql 2 = 0(k*- 1 ) , Qj 3 = G(A; e ), 
where the first relation was derived using the asymptotic relation 

r(± + ^ 2 
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which follows from the Stirling formula. These relations show that the left- 



hand-side of ( |3.11| ) behaves in fix as 

-a 4 k 2 + 0{k 1+2 ') 

for k — > oo, and therefore it diverges uniformly in 9 as long as < e < |. 
Consequently, there is K > such that 

special fl fix = special R fix = . 



Let us pass to the relation (fO). Notice first that A — > a 4 as — > oo in 
fix- Furthermore, for < e < § we have 



2 

„2e-l\ L _ /n/'^e-lN „ _ PMlJ2t-\\ n _ /^llAt-l^ 



a = 0(k e ) , ai = C(^ £ - i ) , = (D(k e ~ ) , Cl = 0(k ) , c 2 = 0(F 

Consequently for e < |, the left-hand side of the spectral condition tends to 
a 8 7^ 0, which implies 

spec(Fv) n fix = 

for K large enough. 

The tight necklaces and carpets exhibit a different behaviour. Now we 
replace the intervals J", J'" defined by ( |5.2| ) by J'^ = \k" — i]", k'^ — fj"] and 
J'" = Win ~ Vn'y — Vn] given in a similar way by 



cos(kd)\ < (Ink) e and 



(I 




COS 7T 




V4 





<(lnfc)" e . (5.3) 



It is straightforward to check that 

« ~ 2(7ra)- 1 (ln^)-% r£X ~ 4( 7 ra)- 1 (ln . 

Consider the set fix := [AT, oo) \ U^=i(^n u J'n) with a fixed > 1- If 
k — > oo in this set, the following estimates hold: 

Qq 1 = A lnfc + 0((ln£;) e ) , Q 12 = ^(Ar^ln fc) e ) , Q^ 3 = C((ln/c) e ) , 

with A 7^ 0. These relations show that the left-hand side of ( |3.12| ) diverges 
for e < 1 like (In A;) 2 , uniformly in 9 as k — > oo within fix- By the same 
token, the left-hand side of ( f4.2| ) diverges for e < 1 like (In A;) 4 , uniformly in 
9 1, 92- We infer that there is a > 1 such that 

specifin H fix = specifpv H fix = specifvi H fix = . 
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Now it is easy to estimate the band and gap widths. The points E' n = 
(k' n ) 2 and = (k") 2 around which the bands concentrate are asymptotically 
like c'n 2 and c"n 2 , respectively, by ( |5.1| ). The widths of the excluded intervals 
behave, in the case of a loose connection, as 

\ J n\ ) \ J n\ ) \ J 'n\ ~ COnst ^ • 

Hence the total length B n of the bands contained in the union of the intervals 
J' n , J", and J™ is of order B n ~ const n 1 ^, and the total length L n of the 
adjacent gaps is L n ~ const (n— n 1_e ) « const n. In the case of a tight con- 
nection the band length is estimated instead by B n ~ const n (Inn) _e which 
still gives gap length increasing linearly with n. We sum up our discussion 
with the following result: 

Proposition 5.1 For loosely connected necklaces and carpets the band-to- 
gap ratio satisfies the bound 

— const n 

as n — > oo, with a positive e < | in Examples I and II, and e < \ in 
Example V. On the other hand, for the tightly connected necklaces and carpets 
in Examples III, IV, and VI, we have 

— - ~ const (Inn) -15 

Ln 

as n — > oo ; with any positive e < 1. 
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